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ABSTRACT 



Context. Based on radial velocity measurements it has been found recently that the two giant planets detected around the star 
HD 73526 are in 2: 1 resonance. However, as our numerical integration shows, the derived orbital data for this system result in chaotic 
behavior of the giant planets, which is uncommon among the resonant extrasolar planetary systems. 

Aims. We intend to present regular (non-chaotic) orbital solutions for the giant planets in the system HD 73526 and offer formation 
scenarios based on combining planetary migration and sudden perturbative effects such as planet-planet scattering or rapid dispersal 
of the protoplanetary disk. A comparison with the already studied resonant system HD 128311, exhibiting similar behavior, is also 
done. 

Methods. The new sets of orbital solutions have been derived by the Systemic Console (http://www.oklo.org). The stability 
of these solutions has been investigated by the Relative Lyapunov indicator, while the migration and scattering effects are studied 
by gravitational N-body simulations applying non-conservative forces as well. Additionally, hydrodynamic simulations of embedded 
planets in protoplanetary disks are performed to follow the capture into resonance. 

Results. For the system HD 73526 we demonstrate that the observational radial velocity data are consistent with a coplanar planetary 
system engaged in a stable 2:1 resonance exhibiting apsidal corotation. We have shown that, similarly to the system HD 128311, the 
present dynamical state of HD 73526 could be the result of a mixed evolutionary process melting together planetary migration and a 
perturbative event. 

Key words, planets and satellites: formation, celestial mechanics, hydrodynamics, methods: N-body simulations 



1. Introduction 

Nearly one third of the multiplanet extrasolar planetary systems 
contain pairs of giant planets engulfed in mean motion reso- 
nances. Resonant extrasolar planetary systems have special im- 
portance when studying the formation of extrasolar planetary 
systems as they require a dissipative mechanism operating on the 
planets which is capable of changing the semi-major axis of their 
orbits. Hence, resonant planetary systems may support the plan- 
etary migration scenario, and thus help in answering the ques- 
tion why Jupiter-sized giant planets do not orbit at their forma- 
tion place. According to the classical formation theories such as 
planetesimal accretion followed by core instability or the grav- 
itational instability model for the giant planets, it is generally 
accepted that they formed quite far from their host star (beyond 
the snowline, which is a ~ 3 - 4 AU for a Sun-l ike star), where 
volatile elements can condens e and accrete (e.g. [Hayashi 1981; 
Pollack 1984; Perrvman 2000). On the other hand, observations 
show that a large fraction of Jupiter-sized planets are found very 
close to their host star at distances smaller than a = 1 AU. 

Planetary migration occurs when the planets are still embed- 
ded into a protoplanetary disk and, due to tidal interactions be- 
tween the planets and the disk, their semi-major axis decrease. 
During the migration process, a resonant capture can occur be- 
tween two planets, if certain dynamical conditions are fulfilled. 
Once the capture happens and the orbital decay is sufficiently 
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slow (adiabatic), the resonant configuration is preserved dur- 
ing the migration, and the two planets can travel very close to 
their host star. The very efficiency of this mechanism is clearly 
demonstrated by the well known resonant system around the star 
GJ 876. The giant planets in this system are engaged in a 2:1 
resonance, having very well determined elements due to their 
short orbit al periods (P\ » 30 , P2 ~ 60 days). It is true in gen- 
eral cases (Be auge et al.l 2006|) and also in th e particular case of 
GJ 876 dLee & Pealdl2002t iKlev et al.ll2005l) . that as a result of 
an adiabatic convergent migration process, beside the mean mo- 
tion resonance the orbits of the giant planets also exhibit apsidal 
corotation, (or in other words apsidal resonance), where the os- 
culating orbital ellipses of the giant planets rotate with the same 
mean angular velocity. 

Numerical integrations based on re cently observed o rbital 
data of the planetary system HD 12831 1 dVogt et al.ll2005l) show 
that the giant planets are in a 2: 1 mean motion res onance without 
exhibi ting apsidal corotation. In a recent study, San dor & Klevl 
(2006) offered a mixed evolutionary scenario for this system 
combining an adiabatic migration leading the system into a 2:1 
resonance and a sudden perturbation. This perturbation could be 
either a close encounter between one of the giants and a rela- 
tively small mass (~ 10M e ) planet already existing in the sys- 
tem (which will refer to in the following as planet-planet scatter- 
ing), or the fast termination of the planetary migration due to a 
sudden dispersal of the protoplanetary disk. Both these perturba- 
tions can be strong enough to induce relatively large-amplitude 
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oscillations of the eccentricities and the resonant angles. While 
retaining the mean motion resonance, the planet-planet scatter- 
ing may even break the apsidal corotation between the orbits of 
the giant planets. 

In this paper we will support the above described mixed evo- 
lutionary scenario by the detailed analysis of the n ewly discov- 
ered r esonant system HD 73526. In a recent paper lTinnev et al] 
(2006) reported that the giant planets of HD 73526 are in the 
2:1 resonance. It has been found that only the resonant angle 
©i = 2A2 - A\ - m\ librates, while both ©2 = 2A2 — A\ — w% 
and Azzr = - tj\ circulate. Here, A\ y 2 denote the mean lon- 
gitude of the inner and outer planet, and their periastron 
longitude. However, according to our numerical integ ration, the 
recently published orbital data bv lTinnev et al.l ([2006), as shown 
in Table Q] result in a weakly chaotic and irregular behavior of 
the system. 

Since chaotic behavior is uncommon among the known reso- 
nant extrasolar planetary systems, and may not guarantee the sta- 
bility of the giant planets for the whole lifetime of the system, we 
have searched for regular orbital solutions for the giant planets as 
well. By using the Systemic Console (http : //www . oklo . org), 
we have found new coplanar orbital solutions for the giant plan- 
ets around HD 73526 exhibiting regular behavior. The regular 
nature (and thus the long-term stability) of these orbits have been 
checked carefully by numerical integratio ns and the Relativ e 
Lyapunov indicator chaos detection method ( San dor et al.l2 004). 
Numerical integrations based on our newly derived orbital data 
exhibit very similar behavior of the eccentricities of the planets 
in HD 73526 to those around HD 1283 1 1 . 

This paper aims at a complete analysis of the resonant system 
HD 73526 including the stability investigation of the new orbital 
solutions and modeling the formation of the resonant system 
HD 73526. The paper is structured as follows: first we present 
four sets of new orbital data. Then we investigate the stability of 
these orbital solutions by using the Relative Lyapunov indicator 
chaos detection method, and finally, we offer possible evolution- 
ary scenarios for the system, which are similar to the studied 
ones in the case ofHD 128311. 



2. Original orbital data and their stability 

It has been m entioned also in the case of HD 128311 by 
IVogt et all d2005l) . that a kinematic orbit fit for orbital parame- 
ters in the case of a strongly interacting resonant pair of giant 
planets based only on two unperturbed ellipses, often results in 
a destruction of the system within a few thousand years. This is 
true for HP 73526 as well. 

Thus lTinnev et alj d2006) have also carried out a self consis- 
tent three-body dynamical orbit fit to the observed velocities (as 
listed in Table 1). They found this dynamical fit to be stable over 
1 Myr integration time. Based on these results of iTinnev et alj 
(2006), we have also performed numerical integration of the sys- 
tem. We have found that the giant planets around HD 73526 are 
indeed in a 2:1 mean motion resonance having libration in the 
resonance variable ©1 with an amplitude ~ 90°. 

However, the behavior of the eccentricities is irregular, see 
Figure 1. This means that the system is chaotic, which of 
course may not exclude a practical stability of the giant planets. 
However, based on the behavior of the observed resonant sys- 
tems to date, we would prefer orbital solutions exhibiting regular 
orbits showing a clear dynamical stability for the whole life-time 
of the system. 

By using th e Relative Lyapunov indicators (RLI, 
Sandor et al. 2004), we have mapped out the stability properties 
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Fig. 1. Irregular behavior of the eccentricities of the giant plan- 
ets in the system HD 73526 using as initial conditions the data 
provided by the original dynamical orbit fit (Table [T]). The upper 
(lower) curve shows the behavior of the eccentricity of the inner 
(outer) planet. 

Table 1. Original dynamical orbit fit of HD 73526 



planet 


mass [Mj] 


a [AU] 


e 


M [deg] 


tzr [deg] 


inner 


2.9 


0.66 


0.19 


86 


203 


outer 


2.5 


1.05 


0.14 


82 


13 



rms [m/s] offset velocity [m/s] 



1.57 



7.9 



-29.96 



of the phase space in the close neigh borhood of the o riginal 
orbital data set, in a similar way done bv lBois et alj ([2003) in the 
case of HD 160691. We have calculated the stability properties 
of the a\ - d2, e\ - e2, M\ - M2, and vj\ - W2 parameter planes, 
where a is the semi-major axis, e is the eccentricity, M is the 
mean anomaly and m is the longitude of periastron of one of the 
giant planets, while the indices ' 1 ' and '2' refer to the inner and 
the outer planet, respectively. 

In Figure 2 we display the stability structures of the parame- 
ter planes for the semi-major axis and the eccentricities. During 
the calculation of a particular parameter plane the other orbital 
data have been kept fixed. This implies that when investigating 
the parameter plane of the semi-major axes a\ -02, the other or- 
bital data such as the eccentricities (ei,<?2), the mean anomalies 
(M\,M2), and the arguments of the periastrons (cti,^) have 
been kept fixed to their original values (Table 1). On each pa- 
rameter plane we have marked the stable regions by white, the 
weakly chaotic regions by grey and the strongly chaotic regions 
by black colors. The actual values of the corresponding orbital 
data are also shown by a symbol '+' on each parameter plane. By 
studying carefully the different stability maps (also the not dis- 
played but calculated cases of the parameter planes tu\ - W2 and 
Mi -M 2), we conclude that the orbital data given bv lTinnev et al.l 
(2006) are embedded in a weakly chaotic region. We stress that 
this does not imply automatically the instability of their fit, how- 
ever by using these orbital data the system exhibits irregular (or 
chaotic) behavior and may be destabilized at later times. 

We have also studied the stability of the system HD 73526 
with slightly modified initial orbital elements. We have found, 
for instance, that even a small error of the initial semi-major axis 
(6a2 = 0.01) or eccentricity ((5e2 = 0.01) of the outer planet leads 
to the destruction of the system in a few hundred thousand years. 
We note that the above errors (<5«2 and 5e2) are well inside the 
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semi-major axis of the Inner planet [AU] 



Table 2. New dynamical orbit fits of the giant planets around 
HD 73526 



Fir 
rii 


planet mass [Mj] 


n r at n 




m Luegj 


■/TV 1 f\t± Ct 1 

m L ue gl 


l 


inner 


2.42 


0.66 


0.26 


69.8 


206.6 




outer 


2.58 


1.045 


0.16 


163.2 


265.6 


2 


inner 


2.62 


0.66 


0.209 


77.7 


208.5 




outer 


2.56 


1.047 


0.194 


131 


316.2 


3 


inner 


2.415 


0.659 


0.26 


70.7 


202.9 




outer 


2.55 


1.045 


0.107 


170.7 


253.7 


4 


inner 


2.675 


0.66 


0.209 


80.1 


207.1 




outer 


2.53 


1.048 


0.172 


122.4 


327.6 


Table 3. Properties of the dynamical orbit fits for HD 73526. 




Fit No. 


xi 


rms [m/s] 


offset velocity [m/s] 






1 


1.81 


8.36 




-46.67 






2 


1.60 


8.09 




-33.76 






3 


1.87 


8.4 




-49.54 






4 


1.58 


8.04 




-32.33 





0.2 0.4 0.6 

eccentricity of the inner planet 



Fig. 2. Stability maps calculated by using the original dynamical 
orbit fit listed in Table 1 . The structure of the parameter planes 
indicates clearly that the orbital data given by the original fit 
(marked by '+') are in a weakly chaotic region. Here brighter 
areas refer to more stable and darker ones to unstable regions. 



error limits of the original orbit fit of iTinnev et al.l d2006l) . being 
Aa 2 = +0.08 and Ae 2 = +0.09. Thus we believe that the use of 
initial orbital data resulting in regular orbits is more convenient 
when investigating the system HD 73526. 

On the other hand, the possibility that the system shows 
chaotic behavior as a result of its formation scenario cannot be 
excluded. It can also happen that a close encounter with a small 
mass body results in a chaotic behavior (characterized by the ir- 
regular variation of the eccentricities for instance), though the 
system itself remains stable for very long time. However, as our 
numerical simulations show, we have not found such a behavior; 
either the system shows regular behavior, or it is destroyed soon 
after the scattering event. 



3. New orbital data and their stability 

In order to find orbital solution for the system HD 73526 ex- 
hibiting regular behavior, we have used the Systemic Console 
(http : //www . oklo . org), and have found some new orbital fits 
for the giant panets. In Table 2 we list four set of them, display- 
ing different behavior. In the case of the first three sets of our 
orbital data, beside the mean motion resonance, the giant planets 
are also in apsidal corotation but with enlarged amplitudes in the 
resonant angles ©2 or Aw, while in the fourth case the apsidal 
corotation is no more present. In all four cases the eccentricities 
of the giant planets show relatively large oscillations. The be- 
havior of the eccentricities is very similar to those found in the 
case of HD 128311 indicating clearly, that the present behavior 
of the system HD 73526 is not likely to be the result of a smooth 
adiabatic migration scenario alone. 

In order to investigate the stability of the newly derived or- 
bital elements (shown in Table 2), we have computed (by using 
the Relative Lyapunov indicator) a series of stability maps, sim- 
ilar to those shown in Figure 2. We have found that in all cases 
the orbital elements are deeply embedded in the stable regions of 
the parameter planes a\ - az, e\ - <?2, vj\ - m^, and M\ - M%. In 
Figure 3 we show the stability structure of the parameter planes 
a\ - fl2 and e\ - ei in the case of Fit 1. The semi-major axes 
and eccentricities corresponding to Fit 1 are marked by the sym- 
bol '+', and they almost lie in the middle of the stability region 
(white regions of the parameter planes). We note that at first sight 
the structure of the parameter planes a.\ — a-i corresponding to 
the original orbital data (shown in Figure 2) and to the orbital 
data of Fit 1 (shown in Figure 3) are quite similar to each other. 
However, in Figure 2 the stable regions are very narrow, the light 
regions (which may indicate stability) are mainly grey and have 
a fuzzy structure indicating clearly the weakly chaotic charac- 
ter of the system. In Figure 3 the light regions are white and they 
have very homogeneous structure showing ordered and therefore 
stable behavior of the system. 

We have also checked by numerical integration whether the 
orbital data shown in Table 2 result in a reliable radial veloc- 
ity curves indicating the motion of the star. Two examples are 
shown in Figure 4, corresponding to the Fit 1 and Fit 4, re- 
spectively. We note that the radial velocity curves of Fit 1 and 
Fit 3 are very similar to each other, while the radial velocity 
curves of Fit 2 and Fit 4 are very similar to the radial velocity 



4 



Zs. Sandor, W. Kley, and P. Klagyivik: Stability and Formation of the Resonant System HD 73526 



100 




0.60 0.65 0.70 0.75 

semi-major axis of the Inner planet [AU] 




0.2 0.4 0.6 

eccentricity of the inner planet 



Fig. 3. Stability maps calculated by using the orbital data pro- 
vided by Fit 1 of Table 2. It can be seen from the parameter 
planes a\ - ai (top) and e\ — e\ (bottom) that the orbital data 
of Fit 1 (marked by '+') are embedded well in the stable white 
region. 
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Fig. 4. Radial velocity curves of the star obtained by numerical 
integration using as initial conditions the orbital data given by 
Fit 1 (top panel) and Fit 4 (bottom banel) of Table 2. It can be 
seen that the calculated curves fit well to the measured radial 
velocity points (black dots). 



[5] show the behavior of the resonant angles ©i and Am, respec- 
tively, corresponding again to Fit 3. 

In the following, when investigating the formation of the sys- 
tem HD 73526, we intend to model a similar behavior provided 
by Fit 1 and Fit 3, in which case the outer planet is a little bit 
more massive than the inner one, and the oscillations in the ec- 
centricities are more moderate than in the cases of Fit 2 and Fit 
4. 



curve given by iTinnev et al.l d2006l) corresponding to their orig- 
inal dynamical orbit fit. In the case of Fit 1 and Fit 3 the mass 
of the inner giant planet is (slightly) smaller than the mass of the 
outer planet, which is also supported by the planetary migration 
scenario. According to full hydrodynamical simulations a cav- 
ity opens between the giant planets, and the inner planet orbits 
in a low density gaseous environment. The outer giant planet 
is still connected to the disk thus it can accrete more material, 
which may resu lt in larger mass for it than for the inner planet 
dKlev et alJl2004T) . 

The similarities between Fit 1 and Fit 3 or Fit 2 and Fit 4 
can also be observed in the behavior of the eccentricities of the 
giant planets. The eccentricities in the cases of Fit 1 and 3 show 
moderate oscillations, while in the cases of Fit 2 and 4 they os- 
cillate with very large amplitudes. One example for the behavior 
of the eccentricities is shown in the top of Figures[5]correspond- 
ing to Fit 3. The other two (middle and bottom) panels of Figure 



4. Formation of the system HD 73526 by smooth 
inward migration 

Since the present behavior of the planetary system around 
HD 73526 is very similar to the behavior of the system 
HD 128311, we shall investigate in this section the possibility 
whether this behavior could be the result of a mixed formation 
scenario combining a smoot h migration with a sud den pertur- 
bative effect as suggested by San dor & Klevl d2006l) in the case 
of HD 128311. We first model the evolution of the system by 
applying a slow inward migration. Having estimated the basic 
characteristics of the planetary migration, we intend to study the 
formation of the system taking into account (i) a fast dispersal 
of the protoplanetary disk (resulting in the termination of the 
migration of the outer giant planet) and (ii) planet-planet scatter- 
ing. 
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Fig. 5. The behavior of the eccentricities e\ and (top), the res- 
onant angles ©i (middle) and Act (bottom) of the giant planets 
by using the initial orbital data given by Fit 3. 



4.1. A hydrodynamical model for HD 73526 

To study the general feasibility of forming a resonant system 
such as HD 73526 through a migration process of two embed- 
ded planets we performed full hydrodynamic evolutions where 
the accretion disk is treated in a flat two-dimensional approxi- 
mation. The general setup and the applied numerical methods of 
the model s are identical to those use d in the de tailed study of 
GJ 876 by lKley et all (I2005I) . see also|KleyJ (Il999l) . Here we re- 
strict ourselves to a description of the particulars of our model 
forHD 73526, and concentrate of new features. 

For the masses of the two planets we use ni\ = 2.415My 
and ;«2 = 2.55My (see Fit 3 from Table 2), where Mj is 
the mass of Jupiter, and the planets are not allowed to accrete 
any mass during their evolution. The planetary orbits are ini- 



tially circular at distances of r\ = 1 and r2 = 2AU from the 
central star with mass 1.08M G . The flat disk extends from of 
r,nm = 0.2 to r nmx = 4.0AU with an initial surface density profile 
of S(r) = Sor -1 / 2 . Here, the density at ro = 1AU is normal- 
ized such that Eo^/Mo = 7.22 • 10~ 4 which gives a total disk 
mass in the domain of about O.O24M . To make the initial evo- 
lution of the planets more realistic we s uperimpose initial den- 
sity gaps to the unperturbed r "'^-profile dKlevll9 99, 2000). The 
disk is driven by an a-type viscosity with a = 0.01, and the tem- 
perature is specified using a fixed relative vertical thickness of 
H/r = 0.05, no energy equation is solved. The smoothing length 
for the gravitational potential is e = 0.6//, which gives approx- 
imately the correct values of the planetary migration rate when 
compared to the full three-dimensional case. For the torque cut- 
off, i.e. the region around the planet which is excluded in the cal- 



culation of the torques, we use r, 



torq 



0.6R 



Roche' 



In the regular 

(non-resonant) evolution of the planets the exact values of e and 
rtorq should not influence the results too strongly due to the lack 
of material in the extended gap regions surrounding the planets. 
When the eccentricities of the planets increase during the evo- 
lution they periodically move through the disk at their peri- and 
apo-center and these numerical parameter may become more im- 
portant. We have checked that our results do not depended on the 
exact values of e and r torq as long as they are in the above speci- 
fied range. 

The computational domain from r m ,-„ to r max is covered by 
240 radial and 504 azimuthal gridcells where the radial spac- 
ing is logarithmic and the azimuthal equidistant. This gives 
roughly squared gridcells throughout the whole domain. At the 
outer boundary we use damping boundary c onditions where the 
disk p rofiles are relaxed to the initial profile dde Val-Borro et all 
2006). At the inner boundary we use a new type of boundary 
condition where we specify a outflow condition with an az- 
imuthally averaged radial velocity which has the magnitude of 



3 v 

' nan 



(1) 



The typical viscous radial inflow velocity of accretion disks is 
v V isc(r) = 3v/(2r). We use here a 5 times higher value than 
Vvisc(?min) to accomodate an increased clearing of the inner disk 
due to /) the disturbances of the embedded planets and if) the 
vicinity of the central star. Nevertheless, this boundary condition 
will prevent material leaving the inner disk too rapidly which is 
typical ly observed in hy drodynamic evolutions with embedded 
planet dCridaet alj|2007b . 

In Fig. [6] we show the distribution of the surface density 
E(r, ip) of the protoplanetary accretion disk after about 1400 
years. The inner planet is located at x — -0.1 \,y = -0.61 and 
the outer one at x — -0.48, y = -1.37. In contrast to previous 
models where the two pl anets orbit in an inner cavity of the disk 
without any inner disk dKlev et alJ l2004). here the inner disk has 
not been cleared and is still present due to the more realistic in- 
ner boundary condition given in Eq. (Q]). 

As a consequence the inner planet still is in contact with the 
inner disk which exerts gravitational torques on it. These tend 
to push the planet outward and will produce a damping of the 
eccentricities. In Fig.|7]we display the evolution of eccentricities 
and semi-major axis of the two planets. Initially, both planets or- 
bit in a wide joint gap where the outer planet experiences only 
the torques of the outer disk and migrates inward, while the in- 
ner disk pushes the inner planet slightly outwards. Capture in 
the 2:1 mean motion resonance occurs at 400 yrs after which 
the planets migrate jointly towards the star. The eccentricities 
increase strongly and reach equilibrium values at approximately 
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Fig. 6. Two dimensional density distribution of a protoplane- 
tary accretion disk with two embedded planets after about 1400 
years. From their initial positions (a\ = 1.0 and «2 = 2.0 AU) 
the planets have migrated inwards and, after capturing into a 2: 1 
resonance, have reached a\ = 0.85 and a-i = 1.36 AU with ec- 
centricities e\ = 0.32 and e2 = 0.08. Within the inner planet 
the inner disk has not been cleared even after such a long evo- 
lution time due to a modified new inner boundary condition for 
the radial inflow velocity. 



1 100 yrs. In contrast additional full hydrodynamical test models 
that do not have an inner disk show a continued eccentricity in- 
crease in particular for the outer planet. We attribute this reduced 
magnitudes of the eccentricities to the damping action of the in- 
ner disk. When the eccentricity of the inner planet increases it 
will periodically have to enter into the inner disk and the gravi- 
tational torques will tend to reduce its eccentricity. Equilibrium 
is then given by the eccentricity driving caused by the inward 
motion of the planets and the damping caused by the disk. After 
resonant capture the resonant angles @i and Am are both librat- 
ing with a libration amplitude of « 35° for ©i and < 10° for 
Am. Hence, this hydrodynamical evolution displays an adiabatic 
migration process which results in apsidal corotation of the two 
planets, where the apsidal lines are always aligned. 

In the long term the eccentricity of the outer planet (ei) re- 
mains at a level of 0.05 to 0.07 and does not decline any fur- 
ther. Additional test computations using different resolutions and 
boundary conditions confirm this trend. The shown hydrody- 
namical simulation represents an explorative study to demon- 
strate that such a damping mechanism of the inner disk may have 
operated. Using this mechanism it might be possible to accomo- 
date a larger radial migration for the reson ant planets than that 
found previously for the case of GJ 876 dKlev et alJ 120051) be- 
cause eccentricties will not be pumped up to very large values. 
This topic will be investigated in more detail in future work. 

4.2. Slow migration modeled by dissipative forces 

Hydrodynamical calculations modeling the migration of plan- 
ets embedded in a gaseous protoplanetary disk require relatively 
large computational time. On the other hand, the damping effect 




Time [Yrs] 

Fig. 7. The evolution of the semi-major axis and eccentricities 
of the two embedded planets for the full hydrodynamical model. 
The inner planet (a l ,e 1 ) is indicated by the darker curves. After 
400 yrs the outer planet captures the inner on in a 2: 1 mean mo- 
tion resonance and they are coupled during their inward migra- 
tion. The eccentricities increase rapidly after resonant capture 
and settle to equilibrium values for larger times. 



of a protoplanetary disk to the orbital evolution of the embedded 
giant planets can also be modeled conveniently in the framework 
of the gravitational N-body problem by using well chosen non- 
conservative forces. These additional forces can be parametrized 
by the migration rate a/a and the eccentricity damping rate e/e, 
or by the corresponding e-folding times r a and r e of the semi- 
major axis and eccentricity of the outer planet, respectively (see 
iLee & Peald 12002; Bea uge et al.l 12006). The relations between 
the damping rates and e-folding times are a/a = —l/r a , and sim- 
ilarly for the eccentricities e/e = -l/r e . We define the ratio be- 
tween the e-f olding times K = T a /T e , (or e/e = -K\a/a\), which 
according to ILee & Peald (120021) . determines the final state of 
the system in the case of a sufficiently slow migration. 

We recall that in the migration scenarios used up to now 
mainly the damping of the outer giant planet has been taken 
into account. As we have mentioned previously, we have also 
supposed that beside the outer disk there is also an inner disk, 
and the planets migrates in a cavity between these disks. The 
presence of the inner disk accelerates the inner giant planet, 
thus forces it to migrate outward and also damps its eccentric- 
ity. This effect can also be modeled by using a (repelling) non- 
conservative force also parametrized by a /a (having positive 
sign for outward migration), and K or e/e. 

In what follows, we present our results obtained in model- 
ing the formation of the system HD 73526 by a slow migration 
process. We have studied three cases. In the first case, beside 
the damping effect of the outer disk, we have taken into account 
the effect of an inner disk. In the second case we have assumed 
that only the outer planet is affected by the outer disk and forced 
to migrate inward. Finally, in the third case we have assumed 
that after the resonant capture the protoplanetary disk dissapears 
gradually inhibiting the further increase of the inner planet's ec- 
centricity. In this latter case only the outer planet has been forced 
to migrate inward. Our aim is to provide an evolution of the sys- 
tem which is in accordance with the eccentricity limits given by 
Fits 1-4. 
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Fig. 8. Evolution of the semi-major axes and the eccentricities 
of the giant planets during different migration scenarios: (i) the 
actions of an outer and an inner disk are modeled, r fl2 = 10 4 
years, r fl] = -5 x 10 4 years and K\ - K2-\^i (Top); (ii) only the 
damping effect of an outer disk is modeled, r fl2 = 10 4 years and 
K = 10, 15, and 20 (Middle); and (iii) the resonant capture takes 
place right before the disk's dispersal obeying a linear reduction 
law between 9 x 10 3 and 1.3 x 10 4 years, r fl2 = 2 x 10 4 years, 
and K=\Q (Bottom). 



In order to compare the results of dissipative N-body cal- 
culations to those obtained by hydrodynamical simulations, at 
the beginning of the migration the giant planets move on cir- 
cular orbits at distances a\ = 1AU and «2 = 2AU around a 
M, = l.O8M mass star. The masses of the giant planets have 
been fixed to m\ = 2.4T5Mj and = 2.55Mj, where Mj is 
the mass of Jupiter. For the inward migration of the outer giant 
planet we have applied T ai — 10 4 years, while for the outward 
migration of the inner giant planet r a , = — 5 X 10 4 years. The 
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Fig. 9. Evolution of the resonant angles ©i and Am during the 
migration scenarios corresponding those displayed in Figure [8] 



ratio between the e-folding times in both cases was K = 10. 
Comparing Figure|7]with the top panel of Figure|8] one can con- 
clude that the presence of an inner disk can be modeled by a 
properly chosen dissipative force. The behavior of the resonant 
angles for this case is shown in the top panel of Figure|9] Initially 
both 0\ and Am deviate slightly from 0° and tend finally to 0°. N- 
body calculations also confirm the results of the hydrodynamical 
simulations that an inner disk can damp efficiently the eccentric- 
ity of an inner planet during a permanent migration. 

Hydrodynamical simulations for HD 73526 also show that 
the lack of the inner disk results in a continuous increase of the 
eccentricity of the inner giant planet. In what follows, we will 
study which value of ratio K is necessary to stop the increase of 
the inner planet's eccentricity when no inner disk is considered. 
(It is known from the hydrodynamical simulations that K lies 
typically between 1 - 10.) 
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We have applied a damping of the semi-major axis of the 
outer planet with an e-folding time T fl2 = 10 4 year using different 
K values (K = 10, 15, and 20). The time evolution of the semi- 
major axes and the eccentricities are shown in the middle panel 
of Figures[8] while the behavior of the resonant angles 9\ and Am 
(with K — 10) is displayed also in the middle panel of Figures [9] 
Studying the behavior of the eccentricities for the different ratios 
K one can see that during the migration process the eccentricity 
of the outer planet is damped sufficiently, while the eccentricity 
of the inner planet is slightly increasing. If the migration is ter- 
minated when the semi-major axes reach their actual values (at 
t « 1.3 x 10 4 year), we need K ~ 15 not to exceed the observed 
eccentricity of the inner planet. (The upper limit for the eccen- 
tricity of the inner planet is around ei max * 0.3, as it has been 
shown by our numerical integrations based on the orbital data of 
Fits 1-4.) On the other hand, from hydrodynamical simulations 
we know that K ~ 1-10, thus although the above values K = 15 
seem to be a bit high, it may be acceptable for thick disks. But 
the K value can be reduced to K\ = K2 — 10 by the presence of 
an inner disk as shown above by the hydrodynamic simulations. 
The evolution of the resonant angles 0; and Am (for K — 10) 
are shown in the middle a panel of Figure [9] After the resonant 
capture, ©1 and Am deviate from 0°, and having reached their 
extrema (~ 20° for ®[ and ~ 50° for Am), they tend to 0° very 
slowly. The reason of this behavior might be that the mass ra- 
tio between the giant planets {rn\jmi - 0.947) lies very close t o 
the critical mass ratio (mi/m^an - 0.95 given by Lee (2004). 
If mi/m.2 > (tni/m2)ctiu the resonant system evolves temporarily 
through an asymmetric apsidal configuration. We note that the 
deviation of the angles ©1 and Am is the largest in the the above 
case when only the outer giant planet is damped by an outer disk. 

It is also clear that the migration of the giant planets should 
terminate when they reach their actual positions. If the termina- 
tion of the migration is a slow process the system reaches a state 
very close to a perio dic solution of the corresponding thr ee-body 
problem studied by Psvc hovos & H adiidemetriou (2005), for in- 
stance. 

In order to avoid the use of the relatively high values for 
K ~ 15 — 20, we have also studied the possibility when the reso- 
nant capture between the giant planets takes place just before the 
dispersal of t he protoplanetary disk. We recall that this scenario 
(proposed by iKlev et aT1l2005l) might have occurred in the case 
of GJ 876, if there is not acting a strong damping mechanism 
on the eccentricities. We have supposed that prior to the reso- 
nant capture the giant planets migrated separately to their orbits 
having a\ = 0.75 AU, and ai = 1.5 AU. Then the outer planet 
migrates inward with an e-folding time of its semi-major axis 
T a — 2 x 10 4 yr with K — 10, and captures the inner planet into 
the 2 : 1 resonance, see the bottom panel of Figure [8] The reso- 
nant capture takes place around t * 5 x 10 3 years, the eccentricity 
of the inner planet grows rapidly, but before exceeding the limit 
<?imax = 0.3, the protoplanetary disk has already disappeared and 
thus the migration is terminated. The disk dispersal happens be- 
tween 9 x 10 3 - 1.3 x 10 4 years, obeying a linear reduction law. 
During this simulation the final value of the inner planet's eccen- 
tricity remained around the critical value e\ < 0.3, however the 
ratio of the e-folding times was only K — 10. 



5. Formation of the system HD 73526 by mixed 
evolutionary scenarios 

In the previous section we have presented three possible migra- 
tion scenarios that can push the giant planets deep into the 2 : 1 
resonance not contradicting to the observed upper limits of the 
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Fig. 10. Top: Behavior of the eccentricities of the giant planets 
during an inward migration followed a sudden dispersal of the 
protoplanetary disk, T aj — 4 x 10 3 and K — 15. The giant plan- 
ets orbit initially at a\ = 1 AU, «2 = 2 AU, and the migration 
is stopped between 6000 and 6050 yrs, applying a linear disk 
dispersal law. The fast termination of the migration results in os- 
cillations of the eccentricities of the planets, very similar to that 
shown in Figure 4. Middle and Bottom: Behavior of the reso- 
nant angles ©i and Am corresponding the above scenario. After 
the fast termination of the migration the resonant angles also be- 
have very similar to the case displayed in Figure 4. 



eccentricities {e\ ~ 0.3, e-i ~ 0.1 ). On the other hand, as our 
dynamic orbit fits (Fit 1-4) in Table [2] show, the eccentricities of 
the giant planets and the resonant angles oscillate with consid- 
erably amplitude. This behavior is similar to that one has been 
observed in the case of the system HD 128311, thus in what 
follows, we intend to model the formation and early evolution 
of the system HD 73526 by using mixed evolutionary scenar- 
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ios melting together inward migration and sudden perturbative 
events mentioned in (San dor & Klevll2"b06l) . 



5.1. Sudden termination of the migration 

Recent Spitzer observations of young stars show that the inner 
part of the protoplanetary disk may be empti ed due to pho- 
toevaporation induced b y the central star (see iD'Alessio et alj 
2005: ICaivet et al.l l2005). Thus when approaching the inner rim 
of such a disk, the inward migration can be stopped rapidly 
dMasset et al.ll2006l) . 

If the inner rim is located at a given radius of the disk (e.g. 
it is stationary) the migration cannot be maintained because in 
this case the the outer disk, which is the driving agent of the 
migration, cannot follow the planet, on its inward drift. 

Based on our experience in modelling the formation of the 
system HD 128311, we have also investigated the above case 
when the termination of the migration (due to the presence of 
a stationary inner rim) happens in a very short timescale, here 
Afdij = 50 years. In this case we apply a faster migration T a = 
4 x 10 3 years and K — 15. The behavior of the eccentricities and 
the resonant angles ©i and Am is shown in Figure[10] It is clear 
that this fast termination of the migration can result in a behavior 
of the eccentricities and the resonant angles, which is similar to 
the observed case shown in Figure 5. 

The reason of this behaviour is that in the beginning of the 
migration process, the resonant angles deviate from 0°, reach 
their extrema, and tend very slowly to 0° (see also the middle 
panel of Figures [9]). If the termination of the migration occurs 
when the resonant angles still deviate from 0°, they begin to os- 
cillate (with an amplitude equal to their values at the end of the 
disk's dispersal) around their equilibrium (which is in this case 
0°). This scenario works also in the cases of slower migration 
rates, however in the case of a faster migration the deviation of 
the resonant angles from 0° is larger, and hence the amplitude of 
the oscillations in e, is noticeably higher. 

We should note that the above fast migration of protoplan- 
ets may not be likely in protoplanetary disks. In thick disks fast 
migration may occur, but in this case is not easy to find relevant 
physical processes leading to the sudden termination of migra- 
tion. However, the outcome of the sudden stop of migration char- 
acterized with the above parameters yields results being in good 
agreement with the behavior of the system by using as initial 
conditions the data of Fit 3. 



5.2. Scattering with a small mass inner planet 

The behavior of the eccentricities of the giant planets around 
HD 73526 is very similar to that observed in the case of 
the system H P 128311 a nd v Andromedae. In the case of v 
Andromedae, iFord et all d2005l) proposed that such a behav- 
ior is most likely the result of a p lanet-planet scat t ering event. 
Studying the system HD 12831 1. Isliidor & Kleyl (120061) have 
demonstrated that the observed behavior of that system may be 
the consequ ence of a planet-plan et scattering as well. In their 
recent paper Tinn ey et alJ ((2006) have also suggested that the 
dynamical behavior of the system around HD 73526, might be 
the result of a dynamical scattering event. In what follows, we 
shall investigate whether the present behavior of HD 73526 can 
be modeled by such an effect. 

We have assumed that in addition to the giant planets, being 
already in a 2 : 1 resonance and migrating towards the hosting 
star, a small mass planet (~ 10M ffi ) is also orbiting close to the 
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Fig. 11. Top: The evolution of the eccentricities of the giant and 
the small mass planets during a planetary migration when the 
inner giant planet captures the small mass inner planet into a 
protective 3 : 1 resonance. This resonant configuration inhibits 
the increase of the small mass inner planet's eccentricity thus 
the system remains stable under the whole migration process. 
Bottom: The behavior of the resonant angles 6\ and Act shows 
that beside the 3 : 1 resonance, the orbits of the inner giant and 
the small mass planet also are in an asymmetric apsidal corota- 
tion. 



hosting star in a quasi circular orbit. (This assumption could be 
realistic as the disco very of a 7.5M© planet around GJ 876 shows 
dRivera et al. 2005)) Initially, the giant planets are far enough 
from the small mass planet thus they do not influence signifi- 
cantly its motion. However, as the giant planets migrate inward 
approaching their present positions, they perturb the motion of 
the small planet. As our numerical experiments show at the cor- 
responding ratio of the semi-major axes, the inner giant planet 
captures the small mass planet into a 3 : 1 resonance. Once the 
capture into the 3 : 1 resonance happens, two different scenarios 
have been detected: (i) the eccentricity of the small mass planet 
grows initially, then oscillates around the eccentricity of the in- 
ner giant planet, (ii) the eccentricity of the small mass planet 
grows during the whole migration process, which can result in 
asurvived possible close encounter between the inner small mass 
and the giant planet. 

In the case (i) the resonant angles of the 3 : 1 resonance 
9\ = 3/ii - A ilm - tb\ - m inn , 82 = 3Ai - A ilm - 2tu\, and 
6*3 = 3/ii - Ai„ n - 2Tn mn librate around ~ 40°, ~ 150°, and 
~ 290°, respectively, while Atzr* = m\ -m mn also librates around 
~ -110° showing that the orbits of the inner giant and of the 
small mass planet are in an asymmetric apsidal corotation. (Here 
the index ' 1 ' refers to the inner giant planet, while the index 'inn' 
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Fig. 12. Top: The evolution of the eccentricities of the giant and 
the small mass planets during a planetary migration when the 
small mass planet is also captured into a 3 : 1 resonance by 
the inner giant planet. In this case the resonant configuration is 
no more protecting, the eccentricity of the small mass planets 
increases during the migration and its orbit can cross the orbit 
of the inner giant planet. Bottom: The behavior of the resonant 
angles 6\ during the migration process. We note that the other 
resonant angles (of the 3 : 1 resonance) are circulating, thus the 
orbits of the inner giant and the small mass inner planet are not 
in apsidal corotation. 



to the small mass inner planet.) This resonant configuration pro- 
tects the small mass inner planet from the close encounter (mak- 
ing impossible the planet-planet scattering), and stabilizes its or- 
bit during the whole migration. The behavior of the eccentricities 
of the planets and the resonant angle 9\ (for the 3 : 1 resonance), 
Aw* are shown in Figure QT| 

In the case (ii) the inner giant planet also captures the small 
mass planet into a 3 : 1 resonance, however, in this case only 
9\ librates around 0° (see the bottom panel of Figure [T2| ). the 
other resonant angles are circulating, thus the orbits of the inner 
giant and the small mass planets are not in apsidal corotation. 
During the migration, the eccentricity of the small mass inner 
planet grows (shown in the top panel of Figure [TZb resulting in a 
more and more elongated orbits. Thus eventually the orbit of the 
small mass planet will cross the orbit of the inner giant planet 
and between the planets a close encounter can happen. After the 
close encounter the small planet is either ejected from the system 
or pushed into a distant orbit from the hosting star. 

Having performed several numerical simulations, we have 
found that the outcome of the resonant capture described in the 
cases (i) and (ii), depends mainly on the initial positions of the 



inner giant and small mass planet, and it is not very sensitive 
on the mass of the small inner planet (ranging between 0.03 and 
Q.Q4Mj llp ). It is very probably that prior to its secular increase, 
e,„„ goes through a sudden change in the family of corotation 
described and analysed re cently also for the 3 : 1 resonance by 
Mich tchenko et al.l (120061) . Thus a detailed analysis would cer- 
tainly be useful to study the mechanism and the probability of 
the resonant capture into the 3 : 1 resonance leading to either 
a protective behavior or a close encounter with the inner giant 
planet and a possible ejection of the small mass planet. However, 
at the present stage of our investigations we restrict ourselves 
only to demonstrate the mechanism that leads to the increase of 
the eccentricity of the small mass inner planet and to a planet- 
planet scattering event. 

In our numerical experiments we varied the mass of the small 
planet between 0.03 and 0.04My„ ; , and the migration timescale 
as well. We have performed several simulations, and have found 
that once scattering happens, in the majority of the cases the 
system remains captured into the 2 : 1 resonance. However, in 
some cases we have observed the complete destruction of the 
resonance, or even escape of the giant planets. We note that if 
the scattering takes place well before the disk dispersal process, 
the large oscillations in the eccentricities may be smoothed out 
due to the damping acting of the disk acting on the eccentric- 
ity of the outer giant planet. If the scattering happens during or 
after the disk dispersal, the oscillations of the eccentricities are 
preserved, and the system reaches its presently observed state. 
We have found that in order to observe the continuous increase 
of the small mass inner planet's eccentricity the migration of the 
giant planets should take sufficiently long time. 

Finally, the question may arise whether it is necessary to 
assume a kind of temporal synchronization between the termi- 
nation of the migration and the scattering event. We have per- 
formed several numerical investigations and found in many cases 
that the inward migrating giant planets represent a perturbation 
strong enough to make the behaviour of the small mass inner 
planet unstable. However, the small planet survived for longer 
times, and the close encounter happened when the migration was 
already terminated. On the other hand there can also be a natu- 
ral temporal synchronization. The migration of the giant planet 
should terminate when they reach their presently observed posi- 
tion. If a small mass planet orbits at a distance from the central 
star, which is needed for the resonant capture, its orbit will be 
captured into a 3 : 1 resonance, and therefore its eccentricity 
pumped up making possible a close encounter with the inner gi- 
ant planet, and finally, its ejection from the system. 

In Figure [T3] which is the magnification of Figure [12] we 
show the typical behavior of the eccentricities of the giant plan- 
ets, and Act during a migration process before and after a planet- 
planet scattering. The giant planets are started initially on distant 
orbits having initial semi-major axes a\ — 2 AU and at — 4 
AU. The e-folding time of the outer giant's semi-major axis is 
t„ 2 = 10 4 years and K = 15. 

The mass of the inner small planet is M,„„ = 0.027 Mj up , 
and it orbits at a- mn = 0.9 AU from the hosting star in a nearly 
circular orbit. The scattering occurs around t sc — 3 x 10 4 years. 
We applied also a linear dispersal process of the protoplanetary 
disk taking place between 2.9 x 10 4 and 3.1 x 10 4 years. Thus the 
oscillations in the eccentricities induced by the scattering event 
are not damped out anymore. Since the giant planets are still in 
a (protective) 2:1 resonance their stability is preserved for the 
whole life-time of the system. We note that in contrary to the 
case HD 128311, during modeling the formation of HD 73526 
we have not observed the breaking of the apsidal corotation. It 
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Fig. 13. Top: The behavior of the eccentricities of the giant 
planets during a planetary migration process before and after a 
planet-planet scattering event taking place around t sc « 3 x 10 4 
years. Bottom: The behavior of Am. The behavior of the eccen- 
tricities and Am is very similar of those of Fit 3. 



is highly probable that due to the large planetary masses and 
relatively small semi-major axes, the system is more fragile than 
HD 128311, and a sufficiently strong perturbation which could 
lead to the breaking of the apsidal corotation may destroy the 
whole system. However, after the scattering event the resonant 
angle Am shows large oscillations around 0°. 

We note that in the above simulations (gravitational three- 
body problem with dissipative forces) the small mass inner 
planet has not been affected by the inner disk. In this case we 
used a model without an inner disk, which may also be realistic 
in some cases. However the damping ratio K = 15 is a bit high; 
smaller K would result in higher eccentricity of the inner giant 
planet. As we have already demonstrated in 4.1, the presence of 
an inner disk can damp effectively the eccentricity of the inner 
giant planet. In order to investigate the effect of the inner disk on 
the small mass planet we have also performed full hydrodynam- 
ical simulations using the same initial parameters as described 
in 4.1, but in addition to the giant planets, a M,„„ = 10M e inner 
planet has also been embedded in the inner disk and started from 
0.4AU. According to the theoretical expectations, its migration 
rate (being in linear regime with no gap opening) is smaller than 
that of a more massive Jupiter-like planet. Thus the small mass 
planet has been captured into a 3 : 1 resonance resulting in an in- 
crease of its eccentricity. In Fig.[14]we display the time evolution 
of the orbital elements of the three planets. At t — 400 orbital pe- 
riods (of the inner giant planet) the two massive planets capture 
each other in a 2 : 1 resonance slowing down the inward mi- 
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Fig. 14. Full hydrodynamical evolution of two embedded mas- 
sive planets and an small additional small mass inner planet. 
Top: The behavior of the semimajor axis of the three planets. 
Bottom: The behavior of the eccentricities. 

gration which nevertheless remains faster than that of the inner 
planet which is consequently captured at a later time. We should 
also note that due to the very long computational time, we have 
not followed the complete evolutionary way of the small mass 
inner planet. On the other hand it is very probable, that the in- 
crease of its eccentricity will result in a scattering with the inner 
giant planet, and finally, ejection from the system. Similarly to 
the eccentricity damping mechanism, the final behavior of the 
inner small mass planet will also be the subject of a more de- 
tailed analysis investigating inner disks with different physical 
properties. However, our results to date clearly show that even 
in the presence of an inner disk the eccentricity excitation mech- 
anism of the small mass inner planet (through resonant capture) 
works well, and enables scattering between the small mass and 
the inner giant planet. 

6. Conclusion 

In this paper we have investigated different evolutionary scenar- 
ios, which may have led to the observed behavior of the reso- 
nant system HD 73526. Accordin g to our numerical i ntegration, 
the original orbit fit provided by iTinnev et alj yOO© results in 
weakly chaotic motion. Thus, we have derived four sets of or- 
bital elements, which fit very well to the observed radial velocity 
measurements, and also provide, through a protective 2:1 reso- 
nance, stable motion for the giant planets during the whole life- 
time of the system. Two of these fits yield very large variations 
of the eccentricities of the giant planets, in one case the giant 
planets are not even in apsidal resonance. In the case of the other 
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two fits the variations of the eccentricities are more moderate, 
and beside the mean motion resonance, the giant planets are in 
apsidal corotation. 

Since the behavior of the giant planet's eccentricities in the 
system HD 128311 is very similar to those in HD 73526, we 
have investigated whether this system may also be the result of 
a mixed formation scenario incorporating a slow migration fol- 
lowed by a perturbative event. Similarly to the formation of the 
system HD 12831 1, we have assumed that the two giant planets 
have been formed far from their hosting star, and migrated in- 
ward due to a planet-disk gravitational interaction. During this 
migration the planets have been locked in a protective 2:1 reso- 
nance. 

Performing full hydrodynamical simulations of two embed- 
ded giant planets in a protoplanetary disk we have found that the 
presence of an inner disk can damp efficiently the eccentricity 
of the inner giant planet resulting in eccentricity values being in 
good agreement with the observations. This mechanism has not 
been investigated yet, and may offer a real possibility to avoid 
the use of very high values of the ratio K, for example in mod- 
eling the resonant system around GJ 876. We have additionally 
found that the effect of an inner disk can be modeled very conve- 
niently by using non-conservative forces. We have investigated 
three convergent migration scenarios that can lead the system to 
the deep 2 : 1 resonance, which are all physically realistic. 

Similarly to HD 128311 we have assumed that the system 
HD 73526 suffered a sudden perturbation near the end of the 
migration of the giant planets. This perturbation could either be 
the quick termination of the migration, which may possibly be 
induced by an inner rim of the disk and an empty region in- 
side of it, as indicated by some observations of young stars, or 
a planet-planet scattering event. We have analyzed the results of 
an encounter of one of the giant planets with a relatively small 
~ 10M e planet, which orbits around the central star at a closer 
distance (a < 1 AU), and which is ejected from the system, or 
thrown to a very distant orbit after the scattering event. 

Investigating the effect of the disk's sudden dispersal, we 
have found that only a very short disk dispersal time Ata sp = 50 
years will induce the observed oscillations in the eccentricities 
and in the resonant angles. The reason of this behavior is that 
during the migration of the giant planets the system tends to 
an asymmetric apsidal corotation with Am + 0°. At moment 
of the sudden termination of the migration Am « 50°, and Am 
will continue to librate with this amplitude around its equlibrium 
(0°). If ©i is also different from 0°, it will also librate around 0° 
after the scattering event. We note however that in order to ob- 
tain the present behavior of the system we need relatively fast 
migration of the giant planets r fl2 = 4 x 10 3 years and K = 15. 

Studying the planet-planet scattering event we have found 
that the inner giant planet always captures the small mass inner 
planet in a 3 : 1 resonance. Then the further evolution of the 
system depends on whether the orbits of the giant and the small 
mass planet are in apsidal corotation after the resonant capture. 
The apsidal corotation protects the orbit of the small mass planet 
from the close encounter with the inner giant planet, thus inhibits 
a planet-planet scattering event. The lack of the apsidal corota- 
tion results in a continuous increase of the small mass planet's 
eccentricity. Thus its orbit crosses the orbit of the giant planet 
and a close encounter can take place between them. Afterwards, 
if the system survives that perturbative event, the eccentricities 
of the giant planets behave qualitatively similar to those of our 
Fit 3, in agreement with the observed radial velocity observa- 
tions of HD 73536. 



We can conclude that based on radial velocity measurements 
to date, it is probable that during the formation and early evo- 
lution of the systems HD 73526 and HD 128311 the adiabatic 
migration, which has brought the systems into a protective 2 : 1 
resonance, has been followed by a sudden perturbative event. 
This perturbative event could be either the fast dispersal of the 
circumstellar disk or a planet-planet scattering. Thus the dynam- 
ics of the system HD 73526 confirm that planetary migration 
followed by sudden perturbations (e.g. planet-planet scattering) 
play an important role in shaping the formation and early dy- 
namics of resonant planetary systems. 
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